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Abstract 

We present simple assumptions on the constraints defining a hard core dynamics for the associ- 
ated reflected stochastic differential equation to have a unique strong solution. Time-reversibility 
is proven for gradient systems with normal reflection, or oblique reflection with a fixed oblicity 
^— I matrix. An application is given concerning the clustering at equilibrium of particles around a 
large attractive sphere. 
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I — I 1 Introduction 

Ph Since the first works of Skorokhod [TT] on existence and uniqueness for strong solutions of re- 
(~| fleeted stochastic differential equations, many authors have investigated this type of equation and 
extended his results on half-spaces to more general domains: convex sets (Tanaka |12]). admissible 
^ sets (Lions-Sznitman [6j), domains satisfying only the Uniform Exterior Sphere and the Uniform 
I— ' Normal Cone conditions (Saisho [7J) or some weaker version of these conditions (Dupuis and Ishii 
I [2]). The question of equilibrium states of the reflected process (construction of time- reversible 
I> initial measures) has also been investigated (see e.g. fW\)- 

All these existence, uniqueness and reversibility results were obtained under some smoothness 
assumptions on the boundary of the domain. Typically, the existence of at least one normal 
inward vector at each point of the boundary is a necessary condition to define the normal reflection 
psj direction. 

In most cases, the domain in which the process has to live is defined by constraints which are 
I physically natural rather than by its geometrical properties as a subset of some Euclidean space. 
For example, consider a system of n identical hard spheres with radius r in R'^. The domain in 
• ^ which they evolve is the set of configurations (xj)i<j<„ satisfying the constraints \xi — Xj\ > 2r 
^ (i.e. the distance between the centers of any two spheres is larger than twice their radius); the 
geometrical description is much more complicated: the complementary set in M of some star- 
convex subset whose boundary can be locally approximated by a tangent sphere and a cone. 

Unfortunately, for reflected processes in dimension larger than three, the geometrical properties 
of the domain are not that obvious from the physical constraints. In the already mentioned nd- 
dimensional example of a finite system of hard spheres, the main part of the paper ^ is devoted 
to proving that the set of allowed ball configurations satisfies the Uniform Exterior Sphere and 
Uniform Interior Cone property. In [S] and [3] too, a meticulous and extensive geometrical study 
has to be performed before the stochastic analysis of the dynamics. 

We present in this note a constraint-based assumption for existence, uniqueness and reversibil- 
ity results for Skorokhod equations. 
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Our aim is to deal with assumptions as simple and physically natural as possible, even if they 
are not the weakest ones. Moreover, we treat the delicate case of oblique reflection. It will be 
the basis of a general approach for analysing stochastic dynamics under constraints. We also 
prove that for well-chosen diffusion coefficients the oblique reflected gradient diffusion admits a 
reversible measure. 

This note is divided in two parts. 

The first part (section]!]) exhibits a new compatibility criterion for constraints. If it is satisfied, 
then the reflected stochastic differential equation admits a unique strong solution, and this solution 
is time-reversible in the case of a gradient system whose reflection direction is consistent with its 
diffusion coefficient. 

In section [3] we present an application. We consider the behaviour of many particles around a 
sphere called planet. The particles are spherical with random radii oscillating between a minimum 
and a maximum value. Their motion is perturbated by collisions into other particles and into the 
planet. The planet also generates a gravitational field which has a smooth attractive influence on 
the particles. We are interested in the equilibrium situation: Are the particles scattered in a wide 
area as a typical configuration, or do they tend to cluster around the planet ? We prove that at 
equilibrium and for low temperature all particles are as close as possible to the planet, all located 
beneath some altitude, with high probability. 



2 Reflected stochastic differential equation under multiple con- 
straints 

We are interested in a process living in the closure of a domain P. This domain is defined by a 
finite set T of smooth M- valued constraint functions on W^: 

P = |x G M'^; /(x) > for each / G -F} . 

T) is an intersection of smooth sets (arbitrary many of them, provided they are in finite number), 
so its boundary is a finite union of smooth boundaries: 

5P = U e P; /(x) = 0} . 

Since we want the process to be reflected on the boundary of P, we have to assume some regularity 
on the functions in T . The reflection at any point x G dT) occurs either in the inward normal 
direction \/f{x) or in an oblique direction depending only on the normal direction. So we have 
to suppose the existence of a direction which is normal to the boundary: V/(x) 7^ for each 
X £ D such that /(x) = 0. We actually assume something more: the flrst derivatives of the 
functions of J- should admit some positive uniform lower bound, their second derivatives should be 
uniformly bounded and, most important, we suppose that the boundary of each single-constraint 
set {x ; /(x) > 0} crosses the boundaries of the other single-constraint sets at "not too sharp 
an angle". To be more precise, we have to exclude infinitely sharp "thorns" which contain a 
point admiting inward normal vectors with opposite directions. This is what we call compatibility 
between the constraints: 

Definition 2.1 Let T be a finite set of M.-valued C"^ -functions on W^. These functions are called 
compatible constraints if 

• P := |x G R'^; /(x) > for each f G -F| is a non-empty connected set ; 

• for each f ^T, inf{|V/(x)|; x G P, /(x) = 0} > and sup{|Z)V(x)|; x G M''} < +00 ; 
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• inf 5(0,Conv(x)) > 

x&dV 



where Conv(x) is the convex hull of the unit normal vectors to the boundaries at point x; 
Conv(x) = 1 g/ |v{[^]| > and ^ cj 

Here and in the sequel, 6 denotes the EucUdean distance in M*^, |y| denotes the Euchdean norm 
of vector y, and \M\ = sup{|My|/|y| ; y E M'^} denotes the norm of the matrix M. Lebesgue 
measure is denoted by dx. 

The next main theorem states that our compatibihty definition provides a convenient assump- 
tion to ensure the existence of a reflected process within a set defined by constraints. In most 
models, for the sake of simplicity, the reflection direction is the inward normal direction on the 
boundary. Here we consider a general oblique reflection, as in the case treated in |4J or in Section 
m We state the result with a fixed deviation 6 *0 from the normal direction. *0 denotes the 
transposed matrix. Normal reflection corresponds to the special case = Id- 

Theorem 2.2 (Existence and uniqueness of the reflected process) 

Let 6 be a fixed d x d invertible matrix and J- be a set of compatible constraints, with T> = 
ri/eJ- ^ /(^) > 0} corresponding subset o/M'^. Ifcr-.V — yR"^^ andh:V — are 
bounded Lipschitz continuous functions on V, then the reflected stochastic differential equation: 

X(t)=x + / cr(X(s))(fW(s) + / b(X(s))(is + V / 6 ^OV f{'X{s))dLf{s) (1) 
^0 Jo /eJ=" ° 

has for each starting point x £ D a unique strong solution in T), where the local times Lf satisfy 

= / I/(X(.))=o dLfis). 
Jo 

In this theorem "Strong uniqueness of the solution" stands for strong uniqueness, in the sense of 
[3] chap. IV def.1.6, of the process X, not of the local times Lf. 

If <T is constant and the drift b is a gradient, then the equation also admits a time-reversible 
measure /i (i.e. the distribution of the solution with initial measure // is invariant under the 
transformation (X(-), (Lj(-))jgjr) — y (X.{T — ■),{Lf{T — ■) — Lf(T))f^jr'^ for each T > 0): 

Theorem 2.3 (Reversibility of the reflected gradient process) 

Let 6 be a fixed dx d invertible matrix and J- a set of compatible constraints. If ^ is a C'^ -function 
on with bounded derivatives, then the solution of 

X(t) = X(0) + 6>W(t) - - /" 6> *6>V$(X(s))ds + V / 6> *6>V/(X(t))dL/(s) (2) 

2 Jo f^^Jo 

admits d/x(x) = lx>{x)e~^^^^dx as a time-reversible measure. 
Remark 2.4: In definition 



2.1 



the condition inf_^ (5(0, Conv) > is equivalent to: 



3/3o>0 VxeaP 3v/0, V/ G ^ s.t. /(x) = v.V/(x) > /3o|v| |V/(x)| 

where the dot denotes the Euclidean scalar product. 

Though this statement is longer and apparently more difficult to obtain than an uniform lower 
bound on the norms of the convex combinations, it is in some sense more intuitive. It states the 
existence of cones (with vertex x, axis v and aperture 2 arccos Po ) which contain all the inward 
normal vectors given by the constraints at point x. The positivity condition ensures that these 



3 



cones do not degenerate into half-spaces. This condition is easier to check in some concrete 
situations (see section [3]). 

Remark 2.5: (Stability of the compatibility property) 

Let be a set of compatible constraints on R'^. 

• If is a d X d invertible matrix, the transformed constraints {f{0 ■); f €z J-} are compatible. 

• If all constraints disregard one of the d coordinates, then induces a set of compatible 
constraints on M*^"^; that is if /(xi, • • • , Xd-i,Xd) = f{xi, • • • , Xd-i,0) for each f m. F and 

each X = (xi, • • • ,Xd) in then |/ : ^ /(x'^) ' ^ ^ compatible. 

The remaining of this section is devoted to the proofs of the above results. We first prove 



remarks 2.4 and 2.5 which will be useful in the other proofs, and then proceed to theorems 2.2 
and [231 



Proof of remark 12.41 

The th ird c ompatibility condition is 3/3o > Vx G &D 5(0, Conv(x)) > /3o- The condition in 

can be rewritten as 3/3o > Vx G &D maxmin \ . — .-. .^'{^ \, - 1 ^ T . /(x) = ol > Bq. 

v^o l|v| |V/(x)| J 



remark 
Thus it 



2.4 



suffices to prove that for each x E dT) 



<5(0,Conv(x)) = maxmin I v.^^; / e .F, /(x) = o| . 

The lower bound on (5(0, Conv(x)) follows from the inequality |y| > y.v > minjgjrj(x)=o |v/(x)| 
which holds for every unit vector v and every y G Conv(x) because families (c/) of non- negative 
numbers summing up to 1 satisfy: 



Since the convex hull Conv(x) is a closed set, it contains an element z with minimal norm: 
|z| = (5(0, Conv(x)). For each / satisfying /(x) = and for each positive e, the convex combination 

jq:^ + e belongs to the convex hull, hence its norm can not be smaller than |z|: 

z V/(x) 

This proves that j— j-.j^y^^ > |z| = (5(0, Conv(x)) and provides the upper bound on (5(0, Conv(x)). 
■ 

Proof of remark 12.51 

Let us prove the compatibility of the set = {g{-) = f{9 •); / G J-"} of transformed constraints. 
Since matrix 6 is invertible, 0~^T> = {y G M'^; V/ G J- f{Oy) > O} is a non-empty connected set 
as continuous image of the non-empty connected set P. 6 also transforms the bounds on the /'s 



into bounds on the (7's. Remark 2.4 with v replaced by 0v provides the existence of some positive 
/3o such that: 

Vx edV 3v ^ 0, V/ G J" s.t. /(x) = v.*0V/(x) > /3o|0v| |V/(x)|. 
Replacing x by Oy we obtain: 

VyG5(0-ip) 3v/0, Vg G -F'' s.t. 5(y) = v.Vg(y) > /?o|0v| \'e-^Vg{y)\ > /3o 'l' '^ff' 



4 



Thanks to remark 2.4 with /3q = |g-i^°|tg| i this proves that is a set of compatible constraints. 

In order to prove the second part of remark [2.5[ we now assume that f{'x.,Xd) = /(x, 0) for 
each f in T and each (x, xj) in M"'. The set P = {x G M"'; /(x) > 0} is equal to P x M where 
P = {z e M'^-^ /(z) > 0} is a non empty connected set as a projection of a non-empty connected 
set. The lower bound on V/ and the upper bound on D^f transfer to / because V/ = (V/, 0) and 
|D^/(x)| = \D'^f{xi, ■ ■ ■ ,Xd-i,0)\- From the compatibility of we also get the existence of a 
positive f3o such that for each x G dV there exists a unit vector v satisfying v.V/(x) > /3o|V/(x)| 
for each function f €z J- vanishing at point x. The last coordinate of V/(x) vanishes, hence 
V = (f 1, • • • , Vd-i) 7^ 0. Since dV = dV x M we obtain the compatibility of the /'s: 

3/3o > Vz G 9P 3v / 0, V/ G s.t. /(z) = v.V/(z) > /3o|v| |V/(z)| 



Proof of theorems [272] and EH] 
The case of normal reflection: we assume here that 6 = Id- According to corollary 3.6 
of [3], equation ([l]) has a unique strong solution as soon as V satisfies the four assumptions 
of the inheritance criterion for Uniform Exterior Sphere and Uniform Normal Cone conditions 
(proposition 3.4 in [3]). We will check these four assumptions in the unusual order (i) (ii) (iv) 
(iii) because some parameter appearing in (iii) depends on a parameter defined in (iv). We use 
the notations: V/ := inf{|V/|(x); x G P, /(x) = 0} and ||L'V||oo := sup{|L»Vl(x); x G W^}. 

Assumption (i): We have to prove that {x G M'^ ; /(x) > O} has boundary in D for each 
constraint /. Let us fix x G P such that /(x) = 0. By definition of the constraint functions, 
V/(x) / 0, that is we can choose an index k such that Va;/(x) ^ 0. For simplicity sake, we assume 
that Vrf/(x) > (the idea easily adapts to k ^ d and to negative partial derivatives). Applying 
the implicit function theorem to the C^-function /, we obtain the existence of a neighborhood 
V of {xi, . . . , Xd^i), a neighborhood U' of Xd and an increasing C^-function h such that the 
C^-diffeomorphism {yi,...,yd) i — > {yi, . . . ,yd-i, f{yi, ■ ■ ■ ,yd)) maps {y G U x U' , /(y) > 0} to 
{{yi,---,yd-i,Zd) (^V xU'; Zd>h{yi,...,yd-i,0)}. Hence, the subset {x G M"*; /(x) > O} has 

boundary in 2?, and its inward normal direction at point x is j^yj-^. 

Assumption (ii): Let us prove that {x G M"'; /(x) > O} satisfies the Uniform Exterior Sphere 
property restricted to P: according to definition 3.1 in [3j, we have to prove that there exists 
some positive a/ such that, for each x G 2? satisfying /(x) = 0, one has 

Vys.t. /(y)>0 (y - x). .'^{j^'j, + — |y - xp > 0. (3) 

Let us fix X G P on which / vanish. Taylor formula gives: 

V/(x).(y - x) + ^(y - x).Z)V(x + c*(y - x))(y - x) = /(y) 
for each y G W^, with some c* G [0; 1] depending on y and x. In particular, for y such that 

/(y) > we obtam (y - x). + 2|v/(x)| - ^^^^^ * ' "-^ ^ \\D^f\\~' 

Assumption (iv): We have to prove the existence of some Pq > such that for each x G dT> there 
exists a unit vector 1^ satisfying l^.V/(x) > /3o|V/| for each constraint such that /(x) = 0. But 



this has already been done in remark 2.4 with /3o = infxeaD d(0, Conv(x)) and Ix = |f| fo^' some 
z with minimal norm in Conv(x). 

Assumption (iii): We have to prove that each set {x G M'^; /(x) > O} satisfies the Uniform Normal 
Cone property restricted to V with constant /3j smaller than /3q/2. Using Taylor formula for the 
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derivative of / yields to V/(y) = V/(x) + Z)^/(x + c*(y — x))(y — x) for some c* G [0; 1] depending 
on y and x. We obtain for x and y on which / vanish: 

V/(x).V/(y) _ |V/(x)| ^ V/(x).D2/(x + c*(y-x))(y-x) 



|V/(x)||V/(y)| |V/(y)| |V/(x)| |V/(y)| 

Since |V/(x)| > |V/(y)| — |-D^/(x + c*(y — x))(y — x)| the right hand side is not smaller than: 

|UV(x + c*(y - x))(y - x)| |V/(x)| |Z)V(x + c*{y- x))(y - x)| 
|V/(y)| |V/(x)||V/(y)| 

fw V/(x) V/(y) ^ , Jl^'/lloo , I , , « ^inir 

that IS \^yjf^ )r|V/( )| — v/ — ~ ■^1' ^ consequence, tor any pf SJO, 1[ one can 

choose a 5j > small enough such that for each x G P satisfying /(x) = and each y £ T> 

v/(y) v/(x) 

|V/(y)nv/(x) 



satisfying /(y) = and |y - x| < 6f one has — — -y.— — -y > Jl- Pj. 



This proves that {x G M*^; /(x) > O} satisfies the Uniform Normal Cone property restricted to V 
with any constant (3f g]0; 1[, in particular with f3f < Pq/2 as requested. 



To complete the proof of theorems 2.2 and |2.3| for 6 = I^, we proceed as in the proof of theorem 
3.3 in [3], replacing the probability measure dfi{x.) = ^lj)(x.)e~^^^^dx in that proof by the {a- 
finite but maybe unbounded) measure // defined by diJ,(x.) = Iii(x)e~*^^^(ix. Girsanov theorem 
yields the density of the distribution of the process with initial measure /i with respect to the 
distribution of reflected Brownian motion with Lebesgue measure as initial measure. Since both 
this density and the distribution of reflected Brownian motion starting from Lebesgue measure 
are time-reversal invariant, we obtain the reversibility of the solution with initial measure ^. 

The case of oblique reflection: Let us check that the results obtained in the normal reflection 
case 6 = Id transfer to the case of any invertible matrix 6. Using the notation = 0~^X, 
existence and uniqueness for equation ([T]) is equivalent to existence and uniqueness for 

x^(t) = x^(o)+ f e~^a{ex.^{s))dw{s)+ f e-^h{ex^{s))ds + y] f 'evf{ex.\s))dLf{s) 

Jo Jo feT^ 

(4) 

in the closure of the set = {y G M'^; V/ G /(^y) > O} with local times satisfying the 



condition Lf{-) = / l/(0x''(s))=o dLf{s). 
Jo 



The transformed coefficients cr^ = 9^a{6 •) and = 0^h{0 •) inherit the boundedness and 
Lipschitz continuity property from a and b. Remark |2.5| provides the compatibility of the set 
of transformed constraints •); / G J-}. Moreover, Q is an equation with normal reflection 

because V{f{0 •)) = ^OV f{6 •). Thus equation Q, and then equation ([l]), have a unique strong 
solution. 

Moreover, if cr = 0, b = -\0^6V<^ and X(0) ~ Ix)(x)e~*Wdx for some C^-function <I> with 



bounded derivatives, then X^ is the solution of equation (4) with cr^ = Id-, = —\ 6\7^{0 •) 
and initial distribution X^(0) ~ lg-i.jy{y)e~^^^^^\det{6)\dy . Thanks to the reversibility result 
obtained for normal reflection, X^ is time-reversible. This implies the time-reversibility of the 
solution of equation ([T]). ■ 



3 An application : cluster of particles near the surface of an 
attractive planet 

As an application of the previous results, let us study the configuration of a large number of 
particles around a fixed sphere called planet. As in [3], we consider spherical hard particles with 
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random radii oscillating between a minimum and a maximum value. Each particle is driven by 
a Brownian motion, and undergoes the influence of the gravitational attraction generated by 
the planet. The motion is perturbated as the particles bump into each other and into the planet. 
What we are interested in is the existence and uniqueness of such a dynamics, and the equilibrium 
distribution of the particles: Do they tend to be scattered in a large area, or are they typically 
close together, all located beneath some altitude, with high probability ? That is, do they form 
some sort of "ocean" around the planet ? Using the results of section [2| we prove in proposition 
|3.3| that the particles eventually tend to cluster at the surface of the planet when the temperature 
(represented by the diffusion coefficient) tends to zero. 

More precisely, the planet is the d-dimensional closed ball -6(0, R) centered at the origin with 
radius R. The space around this planet contains a large number n of particles, which are spheres 
with random positions and random radii. Each one will be represented by the position Xi of its 
center in and the value Xi of its radius. Thus configurations are vectors x = (xi, xi, . . . , Xn, Xn) 
in M'^('^+i). 

To prevent negative radii, we enforce Xi £ [r_, r+] for some fixed values < r_ < r+. 

The random oscillations of the positions of the particles are not on the same scale as the 
random oscillations of their radii. We take this into account via the elasticity coefficient a > of 
their surface. 

We assume that the gravity field f generated by the planet is isotropic: it only depends 
on the norm \x\ (more general cases would be mathematically tractable, but without physical 
meaning). As usual (see e.g. [1]) the gravitational attraction appears through a drift in the 
dynamics: Function is an increasing function which is on ]0; +oo[. The drift decreases with 
the distance, but not too fast (or equivalently the potential energy —ip' increases not too fast) in 
the sense that ip" < and lim infp_>_|_oo Pf'ip) > 0. 

An important example in dimension d = 3 is ^p{p) = C'^*ln(/)) which gives the drift —ip'{p) = 
— corresponding to the usual gravitational acceleration —ip"{p) = 

At temperature > 0, the random motion of particles is modelized by the stochastic differen- 
tial system: 

for i S {1, . . . , n} 

x,{t) = x,(o) + ewi{t) - [ v'{\x,{s)\)^^{s) ds 







Jo R + Xi fr^Jo 



Xi — Xj 
Xi + Xj 



s) dLij{s) 



X,{t) = 1,(0) + eam{t) - ^^Lfit) - Lfit) + Lr{t) -a^Yl 

i=i 

In this equation, vector (Xj(-), Xj(-))i<j<„ represents the positions and radii of the n particles, the 
Wj's are independent M'^-valued Brownian motions and the Wi^s are independent one-dimensional 
Brownian motions, also independent from the Wj's. The amplitude of the Brownian oscillation of 
the positions depends on temperature 9, while the amplitude of the radii oscillation depends on 
both the temperature 6 and the surface elasticity a. Since ip' is positive, the drift of Xi is always 
directed toward the origin, i.e. toward the planet as expected. The local time represents the 
repulsion received by the i*^ particle when it collides with the planet (impulsion away from the 
origin, in the direction of the unit vector "^v ), and the local times L,,- represent the collisions 
between particles, which tend to move the involved particles away from each other: unit direction 
— :J-. Collisions between particles are symmetric (Ljj = Ljj). These local times also appear 

in the dynamics of the radii, because particles are like bubbles, and collisions have a defiating 
effect on them (the radii decrease). The local times Lf and L~ are here to comply with the 
condition Xi G [r_ , r_|_] and give a positive (resp. negative) impulsion to the radius if it becomes 
too small (resp. too large). The impulsions are only given on the boundary of the set of allowed 
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configurations, therefore the if 's, i^'s, 's and Ljj's should satisfy: 



(4) 



for i,j G {1, . . . ,n} 

Lf{t)- 



^\X,{s)\=R+Xi(s) '^^fi^) 



Lfit) 



dLfis) 



L-it) 



The set of constraints is therefore: 

• //^(x) = — (i? + Xi)"^ > for 1 < i < n (the particles do not intersect the planet); 

• /j^(x) = r-i- — Xj > for 1 < i < n (radii are smaller than the maximum value); 

• fj^ (x) = Xj — r_ > for 1 < i < n (radii are larger than the minimum value) ; 

• fij{'x.) = \xi — Xjp — {xi + XjY > for i / j in {1, 2, . . . , n} (particles do not overlap). 
Proposition 3.1 

{/i^) ft^ /i^i 1 < ^ < f^} U {fij', 1 < ^ < i < ?^} is a set of compatible constraints on W^'^'^^'^\ 

Let V = nr=i [ifl'rHRX) n {ftr'm) n (/+)-Hffi;) n n,y.(^-)-nK;))- We want to study 
the effect of the attraction of the planet on a "typical" particle configurations, i.e. at equilibrium. 

Proposition 3.2 

If if is an increasing C"^ -function on ]0;+oo[ satisfying ip" < and lim'mf pip' (p) > 0, then 
equation (£g,£g) has a unique strong solution, which is a V-valued process. 

1 V^?T. /I I \ 

The measure Ix)(x)e~s^ ^'^^ "^^'^'''(ix is a time-reversible measure for the solution. For 9 small 
enough, this measure is finite thus the solution admits a time-reversible probability measure: 



pe{dx) 



dy 



-l2?(x)e~5^S-i^(l^'l)(ix 



Our aim, once the existence and uniqueness of the dynamics is proved, is to verify that, at low 
temperature, all particles cluster around the planet with high probability. In other words, there 
exists with high probability an interface between two regions around the planet: no particle over 
a certain altitude, and beneath this altitude a particle density so high that one cannot add any 
particle more (see figure [l]) . 




Figure 1: A configuration with an interface between high particle density and empty space. 
Proposition 3.3 

For each positive e, let be the set of configurations which do not pack into a minimal volume: 

= {x € 3y €z D 3k < n s.t. \/i ^ k yi = Xi and \yk\ < \xk\ — s} 
The probability that occurs at equilibrium tends to zero as the temperature tends to zero: 



lim pe{Ae) = Q. 
e-s>o 
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The end of the paper is devoted to the proofs of the three above propositions. 

Proof of proposition |3.1| 

The constraints in J-" = {fij ; 1 < z < j < n} U {/j^, //^ ; 1 < i < n} are and the cor- 
responding set D of possible configurations is obviously a non-empty connected set. The first 
derivative of each constraint function is uniformly positive on its vanishing set because: 

. V/f(x) = 2(0,...,0,2;i,-(i? + x,),0,...,0) 

if /«(x) = i.e. \xi\= R + Xi then |V//^(x)| = 2V2{R + Xi) > 2y/2{R + r_) > ; 

• V/+(x) = -V/r(x) = -(0,..., 0,1,0,..., 0) + d2 + 1) - l)th coordinate) ; 

• V/ij(x) = 2 (0, . . . ,0,Xi - Xj, -{xi + Xj),0, .. .,{),Xj - Xi, -{xi + Xj),0, ... ,0) 
if fij{x) = i.e. \xi — Xj\ = Xi + Xj then |V/jj(x)| = 4(xj + Xj) > 8r_ > 0. 

In order to check the condition inf d(0, Conv(x)) > 0, in the form given in remark 2.4l we have 
to find some positive /3o, and some non- vanishing vector v depending on x € &D, such that 

V/ G ^ s.t. /(x) = v.V/(x) > /3o|v| |V/(x)| 

From an intuitive point of view, v is the "shortest way to go back" into V from the point x on 
the boundary of T>, i.e. the quickest way for colliding particles to go apart, for particles with 
maximum (resp. minimum) radius to become smaller (resp. larger) and for particles touching the 
planet to go away. will denote the indices of these globules: = {i s.t. = i? + Xj}. 

Intuitively, the best way to separate colliding particles is to move them away from the center of 
gravity of the cluster, that is to give each center Xi an impulsion in the direction Xi — SjeCi 
where Cj C {1, . . . ,n} is the cluster of colliding particles around Xi (i.e. Q is the set containing 
i and all indices connected to i in the graph constructed on the vertices {1, . . . ,n} by the edges 
j ~ j' <;=^ \xj — Xj'\ = Xj + Xji). Similarly, the best way for particles touching the planet to 
go away is for each center Xi with i G to receive a small impulsion proportional to Xi (this 
impulsion will also separate clusters of colliding particles). So a convenient v should be: 

Xi--^^x, ifC,nC« = ^ f '^-/2 ifxi = r_ 

j^Q^ and Vi = < — r_/2 if Xi = r+ 

Xi ' if n 7^ i otherwise 

Let us prove that the above vector v satisfies the desired inequalities: 

■n I p^- .1. V. V/f(x) R + x^ V R 

• II \Xi\ = R + Xi then Vi = Xi hence: v.- p , , , = — — > —= 

' ' |V//^(x)| V2 ^/2 - V2 

■f- V/^W . r_ V/r(x) . r_ 

• it Xj = r+ then v.- n— — = —Vi = — ; and it Xi = r_ then v.- — - = Vi = — ; 

|V/+(x)| ^ 2 ' |V/r(x)| 2 

• li \xi — Xj\ = Xi + Xj then Ci = Cj which implies Vi — vj = Xi — Xj thus: 

V/ij(x) _ Xi + Xj _ Vi + Vj ^ r_ 
''■|V/i,-(x)| "4 4-4 

So V. |vj(xii is bounded from below, uniformly in x S dV and f £ such that /(x) = 0. To 



complete the proof of proposition 3.1, it only remains to find a uniform upper bound for |v|: 

" 1 
i=l i; CinC«^0 i; CinC«=0 ^ * iSCi 
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If Ci is any cluster of colliding globules: 

2 

_ |2 



2,u^^2(r^-l)(2tta-l) 



fc=0 



6 



Similarly, if Ci is a cluster with at least one globule at distance R + Xi oi the origin: 

ieCi fc=o 

and the same upper bound holds for a sum over a union of such clusters. Consequently: 

|v|2 <2n(i? + r+)2 + ^r2(n-l)n(2n-l) + ^r2 ^ (^Q _ i)(2yCi - 1) + 

i; CinC«=0 

Since the sum over {i; Q n = 0} is smaller than n{n — l){2n — 1), the norm of v is uniformly 
bounded from above as a function of x. This completes the proof. ■ 



Proof of proposition 3.2 



We use theorem 2.2 with the n{d+\) x n{d+l) diagonal matrix 9 which has n times the sequence 
{6, ... ,9, 9a) as its main diagonal entries. Since the constraints are compatible on ]R"'('^+^), for 
any C^-function # on M"('^+i) with bounded derivatives, the equation: 

X(t) = X(0) + 6>W(t) - - / 6> *6>V$(X(s))ds+ V / 6> *6>V/(X(s))dL/(s) (5) 

2^0 f^^Jo 

has a unique strong solution in the closure of the set D defined by the constraints. Choosing 
^(x) = Y11=i^{\xi\)/9'^ hence V^,$(x) = ^^^ip'{\xi\) and V£^^>(x) = 0, equation becomes 

nt nt pt 

Xdt) ^ X,{0) + eW,{t) - ^'(\Xds)\)T^{s)ds+ 2e^X,{s)dLfn.{s) + y2 2e\X,~-X,){s)dLf^^{s) 

Jo Jo ' i=i'^° 

X,{t)=X,{O) + 0aW,{t) + e^a-^ 2{B. + X,){s)dLfR{s) - L^.+ {t) + L^.-{t) -J^J^ 2(1, + (s) j 

Let us define Lij{-) = 29"^ J^{Xi + Xj){s)dLf^^{s), L+ = 9^a^L^+, Lr = 9^a^L^^ and Lf (•) = 

29"^ Jq{R + Xi){s)dL rB{s). The property Lf{-) = / I/(x(s))=o dLf{s) implies that condition 
(Eg) is satisfied for these new local times. Then the solution of equation (15]) is the solution of 

{£e)- Thanks to theorem 2.3, the measure Ix)(x)e~*(''^(ix = Ix)(x)e~ ^'^^ '^'•'^"'^dx is a time- 
reversible measure for the solution. To complete the proof, let us check that this measure can be 
renormalized as a probability measure for 9 small enough. 

£ - n 

Since £ := liminf^^+oo > 0: Vr/ > 3K > \/p> K ip'(p) > - 

P 

For p > K this integrates into f{p) > f{K) + {£ — r?)(ln p — luK), so for any positive constant c: 



+ 00 



-^{p)pd-^dp < e-=^(^)K^(^-^) / p-^'-'i^+'^-^dp 



K 



+00 



K 



For c large enough to satisfy —c{l — rj) + d <Q the above integral is finite, that is: 

^-c<p{\x\)^^ < +00 

\B(0,R) 
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This leads to the desired normahzation constant, for smah enough to satisfy ^ > c: 



V J {M.<i\B{Q,R)Y 

< e'^'^-^'^l I e-^^(l"l)dx) <+oo 

\jR'i-^\B{0,R) J 
where if = min^ji.^^^ip denotes the infimum on [R; +oo[ of the smooth increasing function ip. ■ 

Proof of proposition |3.3| 

Let ip^ := inf{^"^-^ y G V}. This infimum exists because (p is increasing on ]0; +oo[. We 

fix X € Ag. There exists an aUowed configuration y with all particles at the same position as in x 
except one particle (say, the /c'th) which satisfies \yk\ < \xk\ —£■ Since cp' is a decreasing function: 



.-1 .--I "^l^fcl 



1=1 i=l 

Function ip' admits a limit at infinity: 

n 

• If this limit does not vanish, then > y'da^il) > + £ lim ^'{p) > Vr, '■• 

1=1 

• If lim_|-oo f' = 0, the positivity of i := liminf p'^' {p) implies the existence of a X such that: 

p— >'+00 

\/p>K pip'{p) > — and {R + r+)ip'{p) < - hence {p - R - r+)ip'{p) > - 

Without loss of generality, we can choose K > R + 2nr+ + ne. 

— If, among the j;fc's such that there exists y £ T> satisfying ly^l < \xk\ — £ and yi = Xi for 

n 

i ^ k, at least one has a norm smaller than K, then > ip^ + eip'{K) > ip^ ; 

1=1 

— If no "e-pushable" Xk has a norm smaller than K, then all particles in x are at distance 
at least K from the origin (because K > R + 2nr+ + ne and it is impossible to fill so 
large a volume with only n particles without one of them to be " e-pushable" ) . When the 
Xfc which has the largest norm is shifted at distance R + from the origin, becoming 

n 

yk and giving the configuration y, we have ^2 > + (l^fcl — R — r^)ip' {{xf^D > 

1=1 

£ 



So we obtam > (/?( ) > ip + e for all x G with a positive e equal to eliunp if this limit 

^ — ' — +00 

i=l 

does not vanish and to min{eLp' {K) , |) otherwise. 

n 

An immediate consequence is //^(Ag) < ^6)({x € V; y'd^^tl) > + e'})- 

1=1 

The normalization constant of the probability measure fig is larger than: 

thus /iff (vie) < . 

The denominator does not depend on 9. Dominated convergence theorem ensures that the nu- 
merator converges to zero when 6 tends to 0. So we obtain lim ij,Q{Ag) = 0. ■ 
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